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Abstract: Parts manufacturers use sudden death testing to reduce the testing time of experiments. 
The sudden death testing plan in the literature can only be applied when all observations of failure 
time/parameters are crisp. In practice however, it is noted that not all measurements of continuous 
variables are precise. Therefore, the existing sudden death test plan can be applied if failure data/or 
parameters are imprecise, incomplete, and fuzzy. The classical statistics have the special case of 
neutrosophic statistics when there are no fuzzy observations/parameters. The neutrosophic fuzzy 
statistics can be applied for the testing of manufacturing parts when observations are imprecise, 
incomplete and fuzzy. In this paper, we will design an original neutrosophic fuzzy sudden 
death testing plan for the inspection/testing of the electronic product or parts manufacturing. 
We will assume that the lifetime of the product follows the neutrosophic fuzzy Weibull distribution. 
The neutrosophic fuzzy operating function will be given and used to determine the neutrosophic 
fuzzy plan parameters through a neutrosophic fuzzy optimization problem. The results of the 
proposed neutrosophic fuzzy death testing plan will be implemented with the aid of an example. 


Keywords: neutrosophic fuzzy statistics; fuzzy approach; neutrosophic fuzzy plan parameters; fuzzy 
optimization problem; risks 





1. Introduction 


In life testing experiments, a random sample of items is usually selected, and a single item is 
installed on a single tester for the lot sentencing. This type of testing is costly as the number of testers 
is equal to the number of items that are selected for the testing purpose. Alternatively, the group 
sampling scheme is applied for the testing of more than one item on the single tester. Sudden death 
testing is implemented in groups to reduce the testing/inspection cost of the product. In sudden 
death testing, a random sample of size n is equally distributed to g groups having r items in each 
of the g groups. Reference [1] proposed the sudden test for the Weibull distribution. According to 
Reference [1] “The specimens in each group are tested identically & simultaneously on different testers. 
The 1st group of specimens is run until the 1st failure occurs. At this point, the surviving specimens 
are suspended & removed from testing. An equal set of new specimens numbering is next tested 
until the 1st failure. This process is repeated until one failure is generated from each of the groups”. 
The sampling plan for sudden death testing has been considered by several authors in the literature, 
see for example References [1-3]. 

In the modern era, products are manufactured using advanced technology, which result in high 
quality and reliability. For the testing/inspection of a highly reliable product, it may not be possible to 
wait for the failures of the product for the lot sentencing. The two types of censoring widely applied 
to reduce the testing cost for a highly reliable product, are known as type-I censoring and type-II 
censoring. The testing is said to be type-I censoring if the time of experiment is fixed, and type-II if 
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the number of failures is specified for the testing/inspection of the product. The acceptance sampling 
plans for type-I and type-II censoring when the lifetime follows the Weibull distribution is designed by 
a number of researchers, including, for example, References [4-10]. 

The fuzzy approach is widely used when there is some uncertainty in the proportion of defectives. 
In practice, the experimenter may be indeterminate about the percentage of defectives. In this case, 
the traditional sampling plans can be applied for the inspection of the lot. Hence, a number of 
people have designed efficient sampling plans using the fuzzy approach, including, for example, 
References [11-26]. 

Parts manufacturers use sudden death testing to reduce the testing time of the experiment. 
The sudden death testing plan in the literature can only be applied when all observations of failure 
time/parameters are crisp. According to Reference [27], “all observations and measurements of 
continuous variables are not precise numbers but more or less non-precise. This imprecision is 
different from variability and errors. Therefore, lifetime data are also not precise numbers but more 
or less fuzzy. The best up-to-date mathematical model for this imprecision is so-called non-precise 
numbers”. Therefore, the existing sudden death test plan can be applied if failure data or parameters 
are imprecise, incomplete, and fuzzy. The classical statistics have the special case of the neutrosophic 
statistics when there are no fuzzy observations/parameters. Thus, the neutrosophic fuzzy statistics 
can be applied for the testing of manufacturing parts when observations are imprecise, incomplete, 
and fuzzy. Recently, the authors of Reference [28] introduced the neutrosophic statistics in the area of 
acceptance sampling plan. 

By exploring the literature, and to the best of the author’s knowledge, there is no work on the 
design of neutrosophic fuzzy sudden death testing plan using the Weibull distribution. In this paper, 
we will design an original neutrosophic fuzzy sudden death testing plan for the inspection/testing 
of the electronic product or parts manufacturing. We will assume that the lifetime of the product 
follows the neutrosophic fuzzy Weibull distribution. The neutrosophic fuzzy operating function will 
be given and used to determine the neutrosophic fuzzy plan parameters through a neutrosophic fuzzy 
optimization problem. The results of the proposed neutrosophic fuzzy death testing plan will be 
implemented with the aid of an example. 


2. Design of Proposed Plan 


Suppose that Tni€{TL, Tu} =i =1, 2,3, ..., NN being a random sample from the neutrosophic 
fuzzy Weibull distribution with neutrosophic fuzzy shape parameter my and neutrosophic fuzzy scale 
parameter An. The neutrosophic fuzzy Weibull distribution is defined by: 


Fy(tn;myn, Àn) =1-exp (—(tn/An)”" ), tN >0 (1) 


Some more details on fuzzy based Weibull distribution can be seen in References [29,30]. 

It is assumed that items are tested identically and simultaneously in each group. We propose 
following a neutrosophic fuzzy sudden death testing plan. The quality characteristic beyond the lower 
specification limit L is defined as defective. Thus, the probability of defectiveness is given by: 


PN = Pry{Tw < L} = Fy(L) (2) 
For the given py, the corresponding AnL from Equation (2) is given by: 
WN = —In(1 = PN) = (AnL)™N (3) 


Step 1. Select a random sample nye{n ,ny} and distribute r items into gye{ gL, gu} groups. 
Step 2. Record first failure from ith group (i = 1, 2,...gn) and calculate neutrosophic fuzzy statistic 
On = Din Yini SNE{SL Su}, One{OL, Ou}. 


Symmetry 2018, 10, 403 3 of 10 


Step 3. Accept lot of the product if vy > kyl”, kne{kgz, kau} is a neutrosophic fuzzy 
acceptance number. 


The proposed plan has two neutrosophic fuzzy parameters gye{gz,gu} and kye{kar, kau}. 
The proposed sampling plan is the extension of the plan proposed in Reference [1]. The proposed 
sampling plan reduces to the plan in Reference [1] when g; = gu = g and kg, = kau = k. 
The operational process of the proposed plan is also shown in Figure 1. 







Inspect a random 
sample of size nn 


Distribute r items 
into gn groups 





Accept lot 


Figure 1. Operational procedure of proposed plan. 


By following Reference [1] that Yjx, in vy = iar Yjx follows the i.i.d. neutrosophic exponential 
distribution with neutrosophic parameter Ay “rn, So vy follows the neutrosophic fuzzy gamma 
distribution with neutrosophic parameters (gy,An rn); gne{g_,Zu}- The neutrosophic fuzzy 
operating function (NFOC) will be derived as follows: 


L(pn) = 1 — Gogy (2rnknwN); kye{kaL, kau} (4) 


where Gs, is the neutrosophic fuzzy cumulative distribution function of the neutrosophic fuzzy 
gamma with a neutrosophic degree of freedom ôy. 


Neutrosophic Fuzzy Non-Linear Optimization 


Let 4 and P respectively be the producer’s risk and consumer’s risk. The plan parameters 
gne{gr, gu} and kne{kar, kau} will be determined L(Cagr) > 1-— « at acceptable quality level 
(AQL = pi) and L(CroL) < $ at limiting quality level (LOL = pz). Therefore, the neutrosophic 
fuzzy plan parameters of the proposed plan will be determined by following the neutrosophic fuzzy 
non-linear optimization problem: 


minimize nye{nz,ny } (5) 
L(Cagr) = (1 — Gogn(2rnkwwno)) > 1- a, kne{kt ku}; gne{gL gu} (6) 
L(CroL) = (1 — Gogn (2rnknwni)) < b, kwe{kt ku}; gne{gr, gu} (7) 


The plan parameters gye{¢., gu} and kye{kgz, kau } will be determined through Equations (5)-(7) 
using the grid search method. The following algorithm is used to determine the plan parameters. 
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Step 1. Specify the value of r. 
Step 2. Determine the values of gy and ky using the search grid method through Equations (5)-(7). 
Step 3. Choose the parameters for the plan where indeterminacy interval in gy is minimum. 


The combinations that have smaller values of nye{n,,ny} are selected and reported in Tables 1 
and 2. Table 1 reports for gne{gz,gu} and kye{kgz, kau} when r = 5 and Table 2 reports for 
gnefgL,gu? and kye{kaL, kau} when r = 10. From Tables 1 and 2, we noted the following trend 


in gne{QL, gu} and kne{kat, kau}. 


1. For the fixed values of neutrosophic parameters, gnefgL, gu} and kye{kaL, kau} decrease as r 
increases from 5 to 10. 
2. For the fixed values of neutrosophic parameters, gye{g1, gu} and kye{kaz, kau} decrease as 


LQL increases. 


Table 1. The neutrosophic plan optimal parameters when r = 5. 















































pı p2 SN ky L(Cagr) L(CroL) 
oo Min 19 2473.2 0.952163 0.099984 
Max 21 2717.5 0.962712 0.095018 
eau Min 8 783.6 0.953490 0.099964 
Max 10 953.3 0.975800 0.095052 
WA Min 6 462.9 0.969175 0.099899 
Max 8 592.6 0.988845 0.095073 
oooi ae Min 4 222.1 0.973434 0.099858 
i Max 6 311.3 0.994680 0.095147 
nate Min 3 132.6 0.970167 0.099792 
Max 5 201.2 0.996239 0.095123 
ee Min 3 106.0 0.983213 0.099699 
Max 5 160.8 0.998555 0.095117 
mn Min 2 51.5 0.971999 0.099838 
Max 4 89.4 0.998831 0.095418 
nam Min 2 38.6 0.983592 0.099255 
Max 4 66.9 0.999599 0.095298 
me Min 19 987.8 0.952444 0.099978 
Max 21 1085.3 0.962975 0.095069 
me Min 6 184.6 0.969462 0.099904 
j Max 8 236.3 0.988989 0.095133 
sine Min 4 88.5 0.973691 0.099564 
Max 6 123.9 0.994785 0.095364 

0.0025 
an Min 3 52.7 0.983439 0.099924 
Max 5 80.0 0.996321 0.095083 
hae Min 3 42.1 0.983490 0.099512 
Max 5 63.8 0.998600 0.095342 
ia Min 2 25.6 0.958424 0.099282 
Max 4 44.4 0.997442 0.095051 
TE Min 2 15.2 0.984044 0.099320 
; Max 4 26.3 0.999624 0.096062 
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Table 1. Cont. 
pı p2 SN ky L(CaoL) L(CroL) 
as Min 19 492.7 0.952885 0.099908 
Max 21 541.3 0.963371 0.095048 
Aa Min 8 155.8 0.954335 0.099875 
: Max 10 189.5 0.976377 0.095086 
aes Min 6 91.9 0.969849 0.099547 
ee Max 8 117.5 0.989243 0.095381 
j NA Min 4 43.9 0.974240 0.099688 
; Mar 6 61.5 0.994931 0.095194 
Ae Min 3 26.1 0.971200 0.099658 
i Max 5 39.5 0.996491 0.096118 
PR Min 3 20.8 0.983945 0.099161 
; Mar 5 31.5 0.998668 0.095215 
WA Min 2 74 0.984787 0.099317 
j Mar 4 12.8 0.999658 0.096182 
noi Min 19 245.1 0.953835 0.099928 
; Mar 21 269.2 0.964273 0.095309 
nies Min 8 773 0.955450 0.099925 
Max 10 94.0 0.977123 0.095270 
rr Min 5 39.2 0.950025 0.099569 
at Max 7 52.0 0.982409 0.095945 
ines Min 4 26.1 0.955741 0.099193 
Max 6 36.5 0.988705 0.095461 
Baas Min 3 10.2 0.984642 0.096525 
Max 5 15.3 0.998813 0.096244 
‘aes Min 2 48 0.975190 0.099150 
Max 4 8.3 0.999095 0.096094 
7 Min 2 3.5 0.986232 0.098790 
0.200 Max 4 6.0 0.999729 0.099160 
ne Min 19 80.1 0.957244 0.099165 
Max 21 87.9 0.967418 0.095266 
ne Min 8 25.0 0.959516 0.099143 
Max 10 30.3 0.980098 0.096447 
0.03 aie Min 5 12.6 0.954365 0.096610 
Max 7 16.6 0.985015 0.096118 
: Min 4 8.3 0.960407 0.096094 
0.150 Max 6 11.5 0.990833 0.096298 
at Min 2 2.2 0.954964 0.097352 
; Mar 5 4.5 0.999285 0.098199 
fies Min 18 44.9 0.953778 0.098375 
Max 20 49.4 0.965643 0.096049 
ee Min 8 145 0.963878 0.099439 
Max 10 17.6 0.982585 0.095869 
0.05 7 Min 5 7.2 0.960131 0.097749 
0.200 Max 7 9.5 0.987505 0.096649 
sae) Min 4 47 0.965761 0.095135 
Max 6 6.5 0.992691 0.096047 
PE Min 2 1.2 0.961324 0.080608 
; Mar 6 2.7 0.999915 0.095643 
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Table 2. The neutrosophic plan optimal parameters when r = 10. 
pı p2 8N ky L(Cagr) L(Cror) 
WA Min 19 1236.6 0.952163 0.099984 
: Max 21 1358.7 0.962724 0.095049 
an Min 8 391.8 0.953490 0.099964 
Max 10 476.6 0.975815 0.095115 
à Min 6 231.5 0.969147 0.099791 
: Max 8 296.3 0.988845 0.095073 
0.001 Ti Min 4 111.1 0.973396 0.099670 
; Max 6 155.6 0.994688 0.095301 
See Min 3 66.3 0.970167 0.099792 
Max 5 100.6 0.996239 0.095123 
mre Min 3 53.0 0.983213 0.099699 
Max 5 80.4 0.998555 0.095117 
ae: Min 2 25.8 0.971899 0.099239 
Max 4 44.7 0.998831 0.095418 
R Min 2 19.3 0.983592 0.099255 
: Max 4 33.4 0.999602 0.095903 
nes Min 19 493.9 0.952444 0.099978 
Max 21 542.6 0.963004 0.095147 
st Min 6 92.3 0.969462 0.099904. 
Max 8 118.1 0.989014 0.095364 
aad Min 4 44.3 0.973597 0.099096 
i Max 6 61.9 0.994805 0.095754 

0.0025 

ia Min 3 26.4 0.970450 0.099229 
i Max 5 40.0 0.996321 0.095083 
ro Min 3 21.1 0.983387 0.098644 
i Max 5 31.9 0.998600 0.095342 
Min 2 12.8 0.958424 0.099282 
0.030 Max 4 22.2 0.997442 0.095051 
TA Min 2 7.6 0.984044 0.099320 
: Max 4 13.1 0.999629 0.097616 
Saib Min 19 246.4 0.952809 0.099739 
: Max 21 270.6 0.963430 0.095205 
TE Min 8 77.9 0.954335 0.099875 
: Mar 10 94.7 0.976451 0.095405 
f Min 6 46.0 0.969713 0.099008 
T 0.020 Max 8 58.7 0.989294 0.095848 
f sasi Min 4 22.0 0.974054 0.098745 
; Max 6 30.7 0.994970 0.095980 
iain Min 3 13.1 0.970920 0.098260 
Max 5 19.7 0.996529 0.097257 
joan Min 3 10.4 0.983945 0.099161 
: Max 5 15.7 0.998686 0.096635 
i Min 2 3.7 0.984787 0.099317 
0.100 Max 4 6.4 0.999658 0.096182 
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Table 2. Cont. 


















































pı p2 SN ky L(CaoL) L(CroL) 
Sipi Min 19 122.6 0.953686 0.099588 

i Max 21 134.6 0.964273 0.095309 

MA Min 8 38.7 0.955176 0.099196 

i Max 10 47.0 0.977123 0.095270 

iin Min 5 19.6 0.950025 0.099569 

; Max 7 26.0 0.982409 0.095945 

0.01 Siko Min 4 13.1 0.955231 0.097616 
i Max 6 18.2 0.988843 0.096790 

Aaa Min 3 5.1 0.984642 0.096525 

: Max 5 7.6 0.998847 0.099210 

aie Min 2 2.4 0.975190 0.099150 

Max 3 3.3 0.995249 0.097215 

iai Min 2 1.8 0.985482 0.090371 

i Max 4 3.0 0.999729 0.099160 

— Min 19 40.1 0.956814 0.098132 

Max 21 43.9 0.967745 0.096233 

on Min 8 12.5 0.959516 0.099143 

i Max 10 15.1 0.980490 0.098475 

0.03 sau Min 5 6.3 0.954365 0.096610 
Max 7 8.3 0.985015 0.096118 

N Min 4 4.2 0.958944 0.091311 

: Max 7 6.5 0.995704 0.098411 

a Min 2 11 0.954964 0.097352 

i Max 3 1.5 0.988671 0.098094 

an Min 18 22.5 0.952981 0.096593 

Max 20 24.7 0.965643 0.096049 

sak Min 8 73 0.962650 0.095621 

; Max 10 8.8 0.982585 0.095869 

0.05 AA Min 5 3.6 0.960131 0.097749 
; Max 8 5.3 0.993116 0.097357 

eh Min 4 2.4 0.963476 0.086889 

Max 8 41 0.998501 0.098851 

Eh Min 2 0.6 0.961324 0.080608 

Max 8 1.7 0.999996 0.099397 





3. Application of Proposed Plan 


The application of the proposed plan will be given on the ball bearing guality assurance, as the 
guality characteristic is measurable and as there is a chance that some observations may be fuzzy. 
Suppose that the ball bearing manufacturer is interested in applying the proposed sampling plan, 
butis not certain how to select suitable plan parameters for the testing of his product. Suppose he 
decides to install five items on the single tester. Let AQL = AQL = 0.001, LOL = 0.010, a = 0.05 and 
p =0.10. From Table 1, we have gne{3, 5} and kye{106.0, 160.8}. Hence, he can select g between 3 
and 5. Suppose experimenter decides to select a random sample size n = 25 and g = 5. 


Step 1. Select a random sample 25 and distribute five items into five groups. 

Step 2. Perform sudden death testing and note down the first failure from each of the five groups 
(i=1, 2,...,5). The number of first failures from the five groups are Yı = [220,230], 
Y2 = [300, 320], Y3 = 285, Y4 = [155, 165] and Y5 = [365, 375]. The lifetime of ball bearing follows 
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the neutrosophic Weibull distribution with parameter mye{2, 2} and lower specification limit 
L = 200. The statistic vy is calculated as: 


8N 
un = )_ YIN; gne{3, 5}, vne{106.0, 160.8} 
i=1 


vn = [2202, 2302] + [3002, 3202] + [2852, 2852] + 11552, 1652] + 13652, 3752] = [376875, 404375]. Now 
ky LN €{4240000, 6432000}. As vy < ky LN, so reject the lot of ball bearing product. 


4. Comparison Study 


In this section, we will compare the efficiency of the proposed plan using the neutrosophic 
statistics and the sudden death sampling plan proposed in Reference [2] using the classical (crisp) 
statistics in terms of plan parameter g. For fair comparison, we will consider the same plan parameter 
values. By comparing Table 1 of the proposed plan with Table 1 of Reference [2] when r = 5, it can 
be seen that the proposed plan provides smaller values of g as compared to the plan in Reference [2]. 
For easy reference, we present the comparison in gy between the proposed plan and Reference [2] in 
Table 3. For example, when AQL = 0.001 and LOL = 0.002 in Table 3, the proposed plan has a minimum 
value of g = 17 and a maximum value of g = 21, while the plan in [2] provides a crisp value of g = 231. 
Therefore, the proposed plan needs a sample size of n = r*g between 85 and 105. By comparison, the 
existing plan needs a sample size of 1155 for the testing of the same lot of the product. From Table 3, 
we note that the proposed plan has a smaller gy for all combinations of AOL and LOL. Therefore, 
the proposed plan required a smaller sample size for the inspection of the same product. Hence, less 
inspection cost is needed when the proposed sampling plan is implemented. 


Table 3. The comparison of proposed plan with the plan in Reference [2]. 











pi pr Proposed Plan Existing Plan 

8N & 

0.002 [19, 21] 231 

0.003 [8, 10] 154 

0.001 0.004 [6, 8] 115 
0.006 [4, 6] 77 

0.008 [3, 5] 58 

0.100 [18, 20] 34 

0.05 0.150 [8, 10] 16 
0.200 [5, 7] 8 


5. Concluding Remarks 


In this paper, we have designed a neutrosophic fuzzy sudden death testing plan for the 
inspection/testing of the electronic product or parts manufacturing. The NFOC and neutrosophic fuzzy 
non-linear optimization problem are used to evaluate the proposed sampling plan. The proposed plan 
can be applied for the testing of parts when some observations or parameters are fuzzy. The proposed 
plan is the extension of a sudden death plan based on classical statistics. However, the proposed 
plan is more flexible than the plan based on classical statistics. Some tables are given and discussed 
with the help of a ball bearing example. The proposed plan can be applied for testing products in the 
automobile, aerospace and electronics industries. The proposed plan has the limitation that it can only 
be applied when the failure time follows the neutrosophic fuzzy Weibull distribution. The proposed 
plan for some other distribution can be considered as future research. The proposed sampling plan 
using big data can also be considered as future research. 
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